Appljcjtirm transmittal 





A 



UTILITY 




w PATENT APPLICATION 
£ TRANSMITTAL 

>nfy for new non-provisional applications under 37 CFR 1 .53(b)) 



To: Assistant Commissioner for Patents 
Box Patent Application 
Washington D.C. 20231 



Attorney Docket No. Al-Dhahir 1 



First Named Inventor or Application Identifier Naofal Al-Dhahir 



Title Finite-Length Equalization Over Multi-Input Multi- 
Output Channels 



Express Mail Label no EK508640581 US 



- on 0 

u 



APPLICATION ELEMENTS 



Fee Transmittal Form (original and duplicate) 

Specification Total Pages 18 

title 

cross reference to related applications (e.g. provisional application) 

background 

summary 

brief description of the drawings (if filed) 

detailed description 

claims 

abstract 



£3 Drawing(s) 



Total Pages 4 
Total Pages 2 



Ifg Declaration 

yp a. Newiy executed 

0" 1 b. □ Copy from a prior application (37 CFR 1 .63(d)) 

%j (for continuaiions/divisionals with section below filled out) 

H O Deletion of Inventor(s) Signed Statement attached deleting 

fn inventor(s) named in the prior application. 37 CFR 163 (d)(2) 

I 2 *: and 1.33(b). 

Q Incorporation by reference (usable if Declaration is a copy): 

The entire disclosure of the prior application, from which a copy of the oath or declaration 

■s: is supplied, is considered as being part of the disclosure of the accompanying application 

«sj is hereby incorporated by reference herein. 



Q Other 



ACCOMPANYING APPLICATION PARTS 



[X| Assignment 

|^| Recordation form 

^ Power of Attorney 

K| Postcard 

Q Small entity statement 

l~l Certified copy of priority documents 

^ Information disclosure statement 

[><] Copies of IDS citations 

□ 37 CFR 3.73(b) Statement 
[El check 

□ Other 



F ^CONTINUING APPLICATION, check appropriate box and supply the requisite information: 

0 Continuation Divisional □ Continuation-in-part (CIP) of prior Application No: 



CORRESPONDENCE ADDRESS 



I | Customer Number or Bar Code Label 



(insert Customer No. or Attach bar code label here) 



Correspondence Address below 



NAME 


Henry T. Brendzel 






ADDRESS 


P.O. Box 574, Springfield, NJ 07081 






COUNTRY 


United States 




FAX (973)467-6589 


SIGNATURE OF APPLICANT ATTORNEY, OR AGENT 


Name 


Henry T. Brendzel 


Reg. No. 26,844 


Telephone 


(973) 467^2025 ^ 




Signature 




Date <o / 








> / 




I hereby certify that this Application is being deposited with the United States Postal Service "Express Mail Post Office to Addressee" service under 
37 CFR 1.10 on the date indicated above and is addressed to the Assistant Commissioner for Patents, Washington D.C. 20231 . 




V t$ Henry Brendzel 




^4 


Date of Deposit (Printed Name of Person Mailing Paper) / 


jJSig nature of Per 


?on Mailin^Paper) 


> 



FEE TRANSMITTAL 

Patent Fees are subject to annual revisions on October 1. 
These are the fees effective November 10, 1998 
Small entity payments must be supported by a small entity statement, 

wioc?, laiyc ciiiity two lilU2>l t/c (JatU. otrt? rUfiTiS r I U'/oO/Uy-TZ. 


Complete if Known 


AoDlication Nnmhpr 


DU71 00/14 


Filina Date 


1 n/8/QQ 
i u/o/yy 


C! ret MamaH Inmntnr 


N a ota I Al-Dhahir 


Examiner Name 




Group/Art Unit 




TOTAL AMOUNT OF PAYMENT 


($) 730 


Attorney Docket ID 


Al-Dhahir 1 



METHOD OF PAYMENT (check one) 



FEE CALCULATION (continued) 



1. The Commissioner is hereby authorized to charge indicated fees 
and other underpayments, and credit overpayments to: 



Deposit Account Number 

500732 

IXI Charge any additional Fee 
Required under 37 CFR 
1.16and1.17 



Deposit Account Name 

Henry T. Brendzel 

fl Charge the issue Fee Set in 37 CFR 
1 .1 8 at the Mailing of the Notice of 
Allowance 



2. Q Payment enclosed: 

□ Check □ Money Order □ Other 



FEE CALCULATION 

FILING FEE 

Fee Description 

Utility Filing Fee 

Design Filing Fee 

Plant Filing Fee 

Reissue Filing Fee 

Provisional Filing Fee 

SUBTOTAL (1) ($) 



Fee Paid 



690 



690 



]?. CLAIMS 



2T°tal 
jblaims 



Independent 
Claims 



Multiply 

Dependent 

Claims 



Claims 
remaining 



20 



Highest 
Paid 



20 



Extra 



□ 



Rate 



18 



78 



260 



SUBTOTAL (2) ($) 



Amount 



3. ADDITIONAL FEES 



Other 
fee 
Other 
fee 



Fee Description 

Surcharge - late filing fee or oath 

Surcharge ~ late provisional filing fee or cover sheet 

Non-English specification 

For filing request for reexamination 

Requesting publication of SIR prior to Examiner action 

Requesting publication of SIR after Examiner action 

Extension for reply within the first month 

Extension for reply within the second month 

Extension for reply within the third month 

Extension for reply within the fourth month 

Extension for reply within the fifth month 

Notice of Appeal 

Filing a bnef in support of an appeal 

Requesting an oral hearing 

Petition to institute a public use proceeding 

Petition to revive - unavoidable 

Petition to revive - unintentional 

Utility issue fee (or reissue) 

Design issue fee 

Plant issue fee 

Petitions to the Commissioner 

Petitions related to provisional applications 

Submission of Information Disclosure Statement 

Recording each patent assignment per property 

(times number of properties) 

Filing a submission after final rejection 

(37 CFR1. 129(a)) 

For each additional invention to be examined 
(37 CFR 1.129(b)) 



SUBTOTAL (2) ($) 



Fee Paid 



40 



40 




Al-Dhahir 2 



Inventor Information 

Inventor One Given Name:: 
Family Name:: 
Postal Address Line One:: 
5 City:: 
State- 
Zip:: 

Citizenship Country:: 

10 

Correspondence Information 

Name Line One:: 
Address Line One:: 
15 City:: 

State or Province:: 
Postal or Zip Code:: 

Application Information 

20 Title Line One:: 

Total Drawing Sheets:: 
Formal Drawings?:: 
Application Type:: 

25 

Continuity Information 

This application is a:: 
> Application One:: 
Filing Date- 
so 

Representative Information 

Attorney Name:: 
Registration Number:: 
Telephone:: 
35 Fax- 
E-Mail:: 



Naofal 

Al-Dhahir 

105 Harter Road 

Morristown 

NJ 

07960 
Iraq 



Henry T. Brendzel 
P.O. Box 574 
Springfield 
NJ 

07081 



Finite-Length Equalization OverMuIti- 
Input Multi-Output Channels 

Yes 



converted utility application of Provisional 

60/158,714 

October 8, 1999 



Henry T. Brendzel 
26,844 

(973) 467-2025 
(973) 467-6589 
brendzel@home.com 



l 
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Finite-Length Equalization Over 
Multi-Input Multi-Output Channels 

Related Application: 

5 This application claims priority from Provisional application No. 60/158,714, 

filed on October 8, 1999. This application is also related to a Provisional 
application No. 60/158,713, also filed on October 8, 1999. 

Background of the Invention: 

io In multi-user communication over linear, dispersive, and noisy channels, the 

received signal is composed of the sum of several transmitted signals corrupted by 
inter-symbol interference, inter-user interference, and noise. Examples include 
TDMA digital cellular systems with multiple transmit/receive antennas, wide-band 
asynchronous CDMA systems, where inter-user interference is also known as 

15 multiple access interference, wide-band transmission over digital subscriber lines 
(DSL) where inter-user interference takes the form of near-end and far-end 
crosstalk between adjacent twisted pairs, and in high-density digital magnetic 
recording where inter-user interference is due to interference from adjacent tracks. 
Multi-user detection techniques for multi-input multi-output (MIMO) systems 

20 have been shown to offer significant performance advantages over single user 
detection techniques that treat inter-user interference as additive colored noise 
and lumps its effects with thermal noise. Recently, it has been shown that the 
presence of inter-symbol interference in these MIMO systems could enhance 
overall system capacity, provided that effective multi-user detection techniques are 

25 employed. 

The optimum maximum likelihood sequence estimation (MLSE) receiver for 
MIMO channels was developed by S. Verdu, "Minimum Probability of Error for 
Asynchronous Gaussian Multiple Access Channels," IEEE Transactions on 
Information Theory, January 1986, pp. 85-96. However, its exponential complexity 
30 increases with the number of users, and channel memory makes its 

implementation costly for multi-user detection on severe-inter-symbol interference 
channels. 
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Two alternative transceiver structures have been recently proposed for 
MIMO dispersive channels as well. These structures, which are widely used in 
practice for single-input single-output dispersive channels, are the Discrete 
Multitone and minimum-mean-square-error decision feedback equalizer (MMSE- 
5 DFE). In the latter category, this includes A. Duel-Hallen "Equalizers for Multiple 
Input Multiple Output Channels and PAN Systems with Cyclostationary Input 
Sequences," IEEE Journal on Selected Areas on Communications, April 1992, pp. 
630-639; A. Duel-Hallen "A Family of Multiuser Decision-Feedback Detectors for 
Asynchronous Code Division Multiple Access Channels," IEEE Transactions on 

io Communications, Feb/Mar/April 1995, pp. 421-434; J. Yang et an "Joint 

Transmitter and Receiver Optimization for Multiple Input Multiple Output Systems 
with Decision Feedback," IEEE Transactions on Information Theory, Sep. 1994, 
pp. 1334-1347; and J. Yang et al "On Joint Transmitter and Receiver Optimization 
for Multiple Input Multiple Output (MIMO) Transmission Systems," IEEE 

15 Transactions on Communications, Dec. 1994, pp. 3221-3231. Alas, the prior art 
does not offer a practical MIMO MMSE-DFE receiver with feedforward and 
feedback FIR filters whose coefficients can be computed in a single computation 
(i.e., non-iteratively) in real-time under various MIMO detection scenarios. 

20 Summary 

An advance in the art is realized with receiver having a multiple number of 
receiving antennas that feed a MIMO feedforward filter that is constructed from 
FIR filters with coefficients that are computed based on environment parameters 
that are designer-specified. Signals that are derived from a multiple-output 
25 feedback filter structure are subtracted from the signals from the multiple outputs 
of the feedforward filter structure, and the resulting difference signals are applied 
to a decision circuit. Given a transmission channel that is modeled as a set of FIR 
filters with memory v , a matrix W is computed for a feedforward filter that results 
in an effective transmission channel B with memory N b , where N b , where B is 

30 optimized so that B op( = argmin e £race(R ee ) subject to selected constraints; R ee 
being the error autocorrelation function. The feedback filter is modeled by 



2 



8:41 AM 9/24/00 



Al-Dhahir 2 

['n °n,xn,w 4 ] ~ B . where n\ is the number of outputs in the feedforward filter, as 

well as the number of outputs in the feedback filter. 

The coefficients of feedforward and the feedback filters, which are sensitive 
to a variety of constraints that can be specified by the designer, are computed by a 
5 processor in a non-iterative manner, only as often as it is expected for the channel 
characteristics to change. 

Brief Description of the Drawing 

FIG. 1 shows the major elements of a receiver in accord with the principles 
10 disclosed herein; 

FIG. 2 presents the structure of elements 23 and 26; 

FIG. 3 is a flowchart describing one method carried out by processor 22; 

and 

FIG. 4 is a flowchart describing another method carried out by processor 

15 22. 

Detailed Description 

FIG. 1 depicts the general case of an arrangement with n\ transmitting 

antennas 11-1, 11-2, ... 11-/7 q , that output signals (e.g., space-time encoded 
20 signals) to a transmission channel, and n 0 receiving antennas 21-1,21-2,... 21-n 0 . 

Each transmitting antenna p outputs a complex-valued signal x p , the signals of the 

n q antennas pass through a noisy transmission channel, and the n 0 receiving 

antennas capture the signals that passed through the transmission channel. The 

received signals are oversampled by a factor of / in element 20 and applied to 
25 feedforward W filters 23. Thus, the sampling clock at the output of element 20 is 

of period T S =T/I, where 7" is the inter-symbol period at the transmitting antennas. 

The transmission channel's characterization is also referenced to T s . 

Filter bank 23 delivers an n-, plurality of output signals {n\ can equal n q for 

example) from which feedback signals are subtracted in circuit 24 and applied to 
30 decision circuits 25 (comprising conventional slicers). The outputs of decision 

circuits 25 are applied to feedback filters 26, which develop the feedback signals. 
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Processor 22 develops the filter coefficients for the filters within elements 23 and 
26 and installs the coefficients in the filters within these elements, as disclosed in 
detail below. 

In the illustrative embodiment disclosed herein, the received signal is 
5 expressed by 

/=1 m=0 

where y[ ;> is the signal at time k at the / h receiving antenna, is the A7? th 
coefficient (tap) in the channel impulse response between the / h transmitting 
antenna and the / h receiving antenna, and n (y) is the noise vector at the / h 
10 receiving antenna. The memory of this path (i.e., the largest value of m for which 
h { n J) is not zero) is denoted by v (,,y) . 

It may be noted that it not unreasonable to assume, that the memory of the 
transmission channel is the same for all ij links (n,xn 0 such links), in which case 

v (,J) =v . Alternatively, the v (,J) limit in equation (1) can be set to that v which 
15 corresponds to maximum length of all of the n, xn 0 channel input responses, 

i.e.,v = max, y v ( ' J) . It may also be noted that all of the variables in equation (1) are 

actually /x1 column vectors, corresponding to the / time samples per symbol in 
the oversampled FIG. 1 arrangement. 

By grouping the received samples from all n Q antennas at symbol time k 
20 into an /i 0 /x1 column vector y k , one can relate y k to the corresponding titX*! 
(column) vector of input samples as follows 

y* = 5X x *-m +, v (2) 

where H m is the MIMO channel coefficients matrix of size njAn { , x k _ m is a size 
A?, x1 input (column) vector, and i\ is a size a? 0 /x1 vector. 

25 Over a block of N f symbol periods, equation (2) can be expressed in matrix 

notation as follows: 
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(3) 



or, more compactly, 

y*+NH:* = ^ X fr+W f -1*-v + ^k+N f -tk • (4) 

The subscripts in equation (4) indicate a range. For example k + N f -V.k means 
5 the range from k + N f - 1 to k , inclusively. 

It is useful to define the following correlation matrices: 

R xy = E\x k+Nf ^. k _ v y k+N ^_ tk \ = R XX H (5) 

Ryy = Elyk+Nf-tkYk+Nf-u] = HR xx H + R ™ * (6) 
R xx = E[ X k+N f -tk-v X k+N f -tk-v] an ^ (7) 
10 R ™ s £[%W r tX + A/ r tt] • (8) 

and it is assumed that these correlation matrices do not change significantly in 
time or, at least, do not change significantly over a time interval that corresponds 
to a TDMA burst (assumed to be much shorter than the channel coherence time), 
which is much longer than the length of the FIR filters in element 23 (in symbol 

15 periods denoted by N f ). Accordingly, a re-computation within processor 22 of the 
above matrices, and the other parameters disclosed herein, leading to the 
computation of the various filter coefficients, need not take place more often than 
once every TDMA burst. Once H, Rxx and R nn are ascertained (through the 
conventional use of training sequences), Rxy and R yy are computed by R XX H* and 

20 HR XX H* + R„„ , respectively. 

In accordance with the principles disclosed herein, element 23 comprises a 
collection of FIR filters that are interconnected as shown in FIG. 2, and the impulse 
response coefficients of element 23 can be expressed by W* = [W 0 * W* -Vt^ ] , 
each having N f matrix taps W,, of size (/ 0 nxn,). That is, W,. has the form: 
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w 



(9) 



w 



wy 



where each entry in W| p qr) is an /x1 vector corresponding to the / output samples 
per symbol. Stated in other words, the matrix W 0 specifies the 0 th tap of the set of 
filters within element 23, the matrix W, specifies the 1 st tap of the set of filters 
5 within element 23, etc. 

Also in accordance with the principles disclosed herein, element 26 
comprises a collection of FIR filters that also are interconnected as shown in 
FIG. 2, and the impulse response coefficients of element 26 is chosen to be equal 
to 



Stated in other words, B 0 specifies the 0 th tap of the set of filters within element 
15 26, the matrix B 1 specifies the 1 st tap of the set of filters within element 26, etc. 

Defining B* = [0 nxnAb B*] , where B* is a matrix of size n, xn,(N f +v) , the 
value of A/ b is related to the decision delay by the equality (A + N b + 1) = (N f +v) . 



10 



[«„, <wJ- B * s [o, - B o) b; ... b;j, (io) 

where B* is expressed by B* == [Bq B* • • B^ ] , with {N b + 1) matrix taps B,. , each of 




The error vector at time k is given by 

^■k = B X k+W,-1:fc-v — W y fc+W/ _i :fc • 

20 Therefore, the n,- x n ( . error auto-correlation matrix is 



(12) 
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R eeS E[E;EJ 

= B'Rj - B*R X/ W - WR^B + W'R/ 

= B*(R Jtx -R xy R^R yx )B + (W*-B*R xy R;)R y/ (W*-B*R xy R;;) 1 ' 
= B*R ± B + G*R ry G 

Using the Orthogonality Principle, which states that E[E k y' k+Ni _ tk ] = 0 , it can 
be shown that the optimum matrix feedforward and feedback filters are related by 

^Opt — BopfRxyRyy 

= B; f R xx H*(HR xx H*+R nn )- 1 (14) 
5 and the n s x n, auto-correlation matrix R ee is 

R ee -£[E*E;] 

= B*R X B (15) 
- B* (R xx - R ^H* (HR XJf H* + R „ ) 1 HR XX )B 
= B-(Ri+H*R^H)- 1 a 

R ee can also be expressed as R ee = B*R" 1 B , where R = R^ +H*R; 1 n H . 

It remains to optimize values for the B matrix and the W matrix such that, in 
response to specified conditions, the trace (or determinant) of R ee is minimized. 
10 The following discloses three approaches to such optimization. 
Scenario 1 

In this scenario, it is chosen to process only previous receiver decisions. 
These decisions relate to different users that concurrently have transmitted 
information that has been captured by antennas 21 and detected in circuit 25. 
is That means that feedback element 26 uses only delayed information and that the 
0 th order coefficients of the filters within element 26 have the value 0. Therefore, 
in light of the definition expressed in equation (10), this scenario imposes the 
constraint of B 0 = I . 

To determine the optimum matrix feedback filter coefficients under this 
20 constraint, the following optimization problem needs to be solved: 
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/77/A7 g R ee = m/A7gB*R" 1 B , subject to B'O = C* 



(16) 



where 



15 



20 



0 = 



n, 

0 



n, 

0 



0 



and C* =r0 nvnA L ] 



It can be shown that the solution to the above is given by 

B 0P ^RCD(0*R0)- 1 C, 
resulting in the error signal 

R ee , m in=c*(OR<D) 1 c. 

If we define the partitioning 

R^ R 12 



R = 



R 12 R 22 



10 where R^ is of size n,.(A + 1)xn,.(A + 1), then 



B opf = 



R11 
R 12 



■n,(A+1) 

R 12 R^ 1 



and 



(17) 



(18) 



(19) 



(20) 



(21) 



(22) 



Ree,mm ~ C R^C, 

where the delay parameter A is adjusted to minimize the trace (or determinant) of 
R eerm „ . Once B opf is known, equation (14) is applied to develop W* pr 

FIG. 3 presents a flowchart for carrying out the method of determining the 
filter coefficients that processor 22 computes pursuant to scenario 1 . Step 100 
develops an estimate of the MIMO channel between the input points and the 
output point of the actual transmission channel. This is accomplished in a 
conventional manner through the use of training sequences. The estimate of the 
MIMO channel can be chosen to be limited to a given memory length, v , or can be 
allowed to include as much memory as necessary to reach a selected estimate 
error level. That, in turn, depends on the environment and is basically equal to the 
delay spread divided by 7 S . 
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10 



15 



20 



Following step 100, step 110 determines the matrices, R nn , Rxx, Rxy, and 
Ryy. The matrix R nn is computed by first computing n = y -Hx and then computing 
the expected value E[n*n] - see equation (8) above. The matrix R,a is computed 
from the known training sequences - see equation (7) above - (or is pre- 
computed and installed in processor 22). In may be noted that for uncorrelated 
inputs, Rxx=l. The matrices R^ and Ryy are computed from the known training 
sequences and the received signal or directly from H and R nn - see equations (5) 
and (6) above. 

Following step 110, step 120 computes R = R~] - H R^H , and the partition 
components, Rn, Ri 2 , and R22, as per equation (20). Following step 120, step 130 
computes R eemm from equation (22) and adjusts A to minimize the trace (or 

determinant) of R eenun , computes B opt from equation (21), and from B opt 

determines the coefficients of the n, x /i, filters of element 26, pursuant to equation 

(10). Step 140 computes W^ f from equation (14), and finally, step 150 installs the 

coefficients developed in step 130 into the filters of element 26 and the coefficients 
developed in step 140 into the filters of element 23. 

A second approach for computing B opt utilizes the block Cholesky 

factorization (which is a technique that is well known in the art): 



R = R xl + H R J, H 



L, 0 




0 




L* 2 


>2 L 3_ 


0 


D 2 _ 


0 


L 3. 



(22) 



LDL\ 



where L 1 is of size n,.(A + 1)x/7,(A + 1). Using the result in equations (18) and (19) 
yields 



& 0 pt - 



'",(A+1) 


r* 

c = 


_L 2 L- 1 _ 


_ 



Ao P t '"' e ^CAopf+ 1 H 



(23) 



and 
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10 



15 



(24) 



R =C*D: 1 C 

= ^(< op(I ---,< (Aop(+1) _ 1 ), 

where the index A opf is chosen (as before) to minimize the trace and determinant 
of R ee .m(« • Usin 9 equation (23), equation (14) can be expressed as follows 

w; p( =B; pf R xx H*(HR x x+R„ n r 1 

= B; p ,(R^+H*R^H)- 1 H*R-: (25) 



d n ,A ol>l e n,A opl 



u n,(A op( +1)-1 e n,(A op( +1)-1 



L" 1 H * R~l 



Yet a third approach for computing B op( and R eemin defines B* = [C* B*] 



and partitions R L into as 
yield 

.[C* B*] 



.^12 ^22. 



where R ^ is of size n, ( A + 1) x n j ( A + 1) , to 



R^ Rf 2 

rx 12 rv 22 

^12 

i\ 12 rv 22 



c 

B 
B 



(26) 



= (R^ - Rf 2 (R 2 X 2 r 1 Rf 2 * ) + (B *+ Rf 2 (R^ 2 r 1 )R^ 2 (B VR^ 2 )- 1 )\ 
where = CR^C and Rf 2 = C*R^ 2 . Therefore, 

B o P f = -Ri2(R22)~ 1 



(27) 



W^=[0^ I, -Rt 2 (R^r 1 ](R, x H*(HR xx H*-HR nn r 1 
R^^R^-R^R^)" 1 ^ (28) 



Scenario 2 

In this scenario it is assumed that users whose signals are received by the 
FIG. 1 receiver are ordered so that lower-indexed users are detected first, and 
current decisions from lower-indexed users are used by higher-indexed users in 
making their decisions, i.e., B 0 is a lower-triangular matrix. The general results of 
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10 



15 



20 



25 



equations (21) and (22) can be applied by setting C* = [Q„ iXnA B*] where B 0 is an 

rifXrij monic lower-triangular matrix whose entries are optimized to minimize 
trace{R e e,min). To this end, a partitioning can be defined where 



R 1 R 2 

R 2 R 3 



(29) 



Rn being the term corresponding to Rn of equation (20), with Ri being of size 
/7 ; Axn ; A , and R 3 being of size n i xn } . Equation (22) simplifies to 

R ee,min = B 0 R 3 B 0 (30) 

It can be shown that, the optimum monic lower-triangular B 0 that minimizes 
frace(R ee ,min) is given by the nomic lower-triangular Cholesky factor of R~ 1 , i.e., 



which yields 



and 



R 3 1 ~ L3D3L3 , 



Qopt _ 1 
D 0 ~ L 3 » 



R ee,min — ^3 



-1 



(31) 
(32) 
(33) 



The result is that B opt = 



'n,(A+1) 



C , as expressed in equation (21), with the 



modified value of R^\ and with 

CX<W B o]- (34) 
A second approach for computing the optimum FIR filter coefficients for the 
FIG. 1 receiver involves computing a standard - rather than a block - Cholesky 
factorization of the matrix R = R'^ + H*R; 1 „H (see the definition following equation 
(15)) in the form LDL*. Then, the coefficients of the element 23 filters is given by 
the nj adjacent columns of L that correspond to a diagonal matrix with the smallest 
trace. Therefore, equations (23) and (25) are used to compute the corresponding 
coefficients, with the understanding that L is now a lower-triangular matrix, rather 
than a block lower-triangular matrix. The equivalence of the two approaches can 
be shown using the nesting property of Cholesky factorization. 
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FIG. 4 presents a flowchart for carrying out the method of determining the 
filter coefficients that processor 22 computes pursuant to scenario 2. Steps 100 
through 120 are the same as in FIG. 3, but the method diverges somewhat in the 
following steps. In step 131 the partition according to equation (20) is developed 
5 for a A that minimizes R ee ,min of equation (33), and control passes to step 141 , 
where B° pf is computed based on equations (31) and (32), followed by a 

computation of B opt based on equations (21) and (34). Following step 141, step 

151 computes W* pf from equation (14), and finally, step 161 installs the 

coefficients developed in step 141 into the filters of element 26 and the coefficients 
10 developed in step 151 into the filters of element 23. 
Scenario 3 

When multistage detectors are employed, current decisions from all other 
users, obtained from a previous detection stage, are available to the user of 
interest. Therefore, suppressing their interfering effects would improve the 
15 performance of the receiver. This detection scenario has the same mathematical 
formulation as scenarios 1 and 2, except that B 0 is now constrained only to be 
monic, i.e., e^e, = 1 for all 0 < / < n,. -1. The general results in equations (21) 

and (22) still apply with C* = [0 nxn A B* ] where B 0 is optimized to minimize 

trace{R ee ,min)- In short, under scenario 3, the following optimization problem is 
20 solved: 

minfrace(B*R 3 B 0 ) subject to e*B 0 e, = 1 for all 0 < / < n, - 1 , (33) 

where R3 is as defined in equation (29). Using Lagranage multiplier techniques, it 
can be shown that the optimum monic B 0 and the corresponding MMSE are given 
by 

25 Bf = R .i & !-' ; 0</<n,.-1. (34) 

e,R- 1 e, 

Thus, the method of determining the filter coefficients that processor 22 computes 
pursuant to scenario 3 is the same as the method depicted in FIG. 4, except that 
the computation of B° 0 pt within step 141 follows the dictates of equation (34). 
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With the above analysis in mind, a design of the filter coefficients of the 
filters within elements 23 and 26 can proceed for any given set of system 
parameters, which includes: 

• MIMO channel memory between the input points and the output point of the 
actual transmission channel, v , 

• The number of pre-filter taps chosen, A/ f , 

• The shortened MIMO memory, /V b , 

• The number of inputs to the transmission channel, n u 

• The number of output derived from the transmission channel, n 0 , 

• The autocorrelation matrix of the inputs, Rxx, 

• The autocorrelation matrix of the noise, R nnj 

• The oversampling used, /, and 

• The decision delay, A . 

It should be understood that a number of aspects of the above disclosure 
are merely illustrative, and that persons skilled in the art may make various 
modifications that, nevertheless, are within the spirit and scope of this invention. 
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Claims: 

1. A receiver responsive to an n 0 plurality of entry points comprising: 
a feedforward filter structure having an n 0 xn s plurality of FIR filters, each 

responsive to a signal that is derived from one of said n 0 entry points and each 
5 developing an output signal that contributes to one of n,- feedforward filter structure 
outputs; 

a feedback filter structure developing n\ feedback signals, the structure 
having an n i xn / plurality of FIR filters, each being responsive to one of rt\ output 
signals; 

10 a subtracter structure that develops n\ signals from signals of said n { 

feedforward filter structure outputs and said n\ feedback signals; and 

decision logic responsive to said m outputs developed by said subtracter 
structure, for developing said n\ output signals. 

15 2. The receiver of claim 1 further comprising a sampling circuit interposed 

between said n 0 plurality of entry points and said feedforward filter structure that 

samples received signal at rate 7 S = y , where / is an integer and T is symbol rate 
of a transmitter whose signals said receiver receives. 

20 3. The receiver of claim 2 where />1 . 

4. The receiver of claim 1 where further including a processor coupled to 
signals applied to said feedforward filter structure, for computing coefficients of 
said FIR filters included in said feedforward filter structure and of said FIR filters 

25 included in said feedback filter structure. 

5. The receiver of claim 4 where coefficients are computed in said 
processor in response to a block of N f symbols. 
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6. The receiver of claim 4 where said coefficients are computed in a non- 
iterative manner. 

7. The receiver of claim 4 where said coefficients are computed with a non- 
5 iterative equation. 

8- The receiver of claim 4 where said coefficients are computed once every 
time interval that is related to rapidity of change in characteristics of transmission 
medium preceding said entry points. 

10 

9- The receiver of claim 8 where said processor installs computed 
coefficients of said FIR filters in said FIR filters following each computation. 

10. The receiver of claim 1 where said FIR filters in said feedforward filter 
is structure form an array of filters that includes one FIR filter connected between 

each of said n 0 input points and said n, outputs. 

11. The receiver of claim 10 where said entry points are antennas and said 
n 0 plurality of antennas receive signals, via a transmission channel, from a 

20 transmitter having a multiple number of transmitting antennas. 

12. The receiver of claim 11 where said transmitter has n,- transmitting 
antennas. 

25 13. The receiver of claim 2 where said plurality of FIR filters in said 

feedforward structure is expressed by matrix W, and W is computed 

by W; pf = B* pf R xy R^ , w; pf = B; p ,R (HR XX H* + RJ' 1 , or 

W op* =B* pf (R; x + H*R^H)" 1 H*R^ , where R xx is an autocorrelation matrix of a block 
of signals transmitted by a plurality of transmitting antennas to said n 0 antennas 
30 via a channel having a transfer characteristic H, R nn is an autocorrelation matrix of 
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noise received by said plurality of n 0 antennas during said block of signals 
transmitted by said transmitting antennas, R xy = R XX H* , = HR x;f H* + R nn , and 

B* p , is a sub-matrix of matrix B* opt , where B opf = argmin e trace{R ee ) subject to a 
selected constraint, R ee being the error autocorrelation function. 

14. The receiver of claim 13 where said plurality of FIR filters in said 
feedback structure is expressed by matrix [l n , 0^ ] - B\ 



15. The receiver of claim 1 where coefficients of the FIR filters in said 
10 feedforward filter are set to results in an effective transmission channel B with 

memory Nb, where N b < v , where B is optimized so that B op/ = argmin B frace(R ee ) 

subject to a selected constraint; R ee being the error autocorrelation function, the 

feedback filter is modeled by [l„ 0 ni xn,N b ] ~ B * > where ri\ is the number of outputs 

in the feedforward filter, as well as the number of outputs in the feedback filter, and 

is the feedforward filter is modeled by W, where W opt = B* p ,R xy R^ , R xy is the cross 

correlation between transmitted signals and signals received by said receiver, and 
Ryy is the autocorrelation of the received signals. 

16. The receiver of claim 15 where said selected constraint is B*fl> = C* , 



20 where O ; 
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17. The receiver of claim 15 where said selected constraint is B*<D = C* , 

and C* =[0 nxnA B* 0 ] , being a monic lower- 
triangular matrix whose entries are optimized to minimize trace(R ee ,min) 

5 18. The receiver of claim 15 where said selected constraint is e,B 0 e,. = 1 , 

where e; is a vector with value 1 in position / and values 0 elsewhere, and where 
B* being a monic lower-triangular matrix whose entries are optimized to minimize 
frace(R e e,m/n). 

io 19. The receiver of claim 13 wherein said plurality of FIR filters in said 

feedback filter structure and in said feedforward filter structure are subjected to 
designer constraints relative to any one or a number of members of the following 
set: transmission channel memory, size of said block, effective memory of the 
combination consisting of said transmission channel; n h n 0l autocorrelation matrix 

is Rxx, autocorrelation matrix R nnj value of factor / in said sampling circuit, and 
decision delay. 

20. The receiver where said matrix W is expressible by 
W = [W 0 W, ■ W WH J , where matrix W g is a matrix that specifies q th tap 
20 coefficients of said FIR filters. 



where O = 



L 0 
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Abstract 

A Ml MO Decision Feedback Equalizer improves operation of a receiver by 
cancelling the spatio-temporal interference effects caused by the MIMO channel 
memory with a set of FIR filters in both the feedforward and the feedback MIMO 
filters. The coefficients of these FIR filters can be fashioned to provide a variety of 
controls by the designer. 
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FIG. 4 
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IN THE UNITED STATES 
PATENT AND TRADEMARK OFFICE 

Declaration and Power of Attorney 

As a below named inventor, I hereby declare that: 

My residence, post office address and citizenship are as stated below next to my name. 

I believe I am an original, first and sole inventor of the subject matter which is claimed 
and for which a patent is sought on the invention entitled Finite-Length Equalization Over 
Multi-Input Multi-Output Channels the specification of which was filed on 10/8/99, as 
provisional application Serial No. 60/158714. 

I hereby state that I have reviewed and understand the contents of the above -identified 
specification, including the claims, as amended by an amendment, if any, specifically referred to 
in this oath or declaration. 

I acknowledge the duty to disclose all information known to me, which is material to 
patentability as defined in Title 37, Code of Federal Regulations, 1.56. 

I hereby claim foreign priority benefits under Title 35, United States Code, 119 of any 
foreign applications) for patent or inventors' certificate listed below and have also identified 
below any foreign application for patent or inventors' certificate having a filing date before that 
of the application on which priority is claimed: 



I hereby claim the benefit under Title 35, United States Code, 120 of any United States 
applications) listed below and, insofar as the subject matter of each of the claims of this 
application is not disclosed in the prior United States application in the manner provided by the 
first paragraph of Title 35, United States Code, 112, we acknowledge the duty to disclose all 
information known to us to be material to patentability as defined in Title 37, Code of Federal 
Regulations, 1.56 which became available between the filing date of the prior application and the 
national or PCT international filing date of this application: 

Provisional Application No. 60/158,714, filed October 8, 1999 

I hereby declare that all statements made herein of my own knowledge are true and that 
all statements made on information and belief are believed to be true; and further that these 
statements were made with the knowledge that willful false statements and the like so made are 
punishable by fine or imprisonment, or both, under Section 1001 of Title 18 of the United States 
Code and that such willful false statements may jeopardize the validity of the application or any 
patent issued thereon. 
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I hereby appoint the following attorneys) with full power of substitution and revocation, 
to prosecute said application, to make alterations and amendments therein, to receive the patent, 
and to transact all business in the Patent and Trademark Office connected therewith: 
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(Reg. No. 27873) 


Thomas A. Restaino 
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Gary H. Monka 
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I also appoint the following as associate attorney(s), with fiiJl power to prosecute said 
application, to make alternations and amendments therein, and to transact all business in the 
Patent and Trademark Office connected therewith: 

Henry T. Brendzel (Reg. No. 26,844 
William Ryan (Reg. No. 26,844) 
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